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Abstract 

We address Newton-type problems of minimal resistance from an op- 
timal control perspective. It is proven that for Newton-type problems 
the Pontryagin maximum principle is a necessary and sufficient condition. 
Solutions are then computed for concrete situations, including the new 
case when the flux of particles is non-parallel. 
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1 Introduction 

In 1686, in his celebrated Principia Mathematical Isaac Newton propounded 
the problem of determining the profile of a body of revolution, moving along its 
axis with constant speed, through some resisting medium, which would minimize 
the total resistance (see [HIE])- Problems of this kind find application in the 
building of high-speed and high-altitude flying vehicles, such as in the design 
of missiles or artificial satellites. Newton has given the correct answer to his 
problem, in the situation of a "rare" medium of perfectly elastic particles with 
constant mass and at equal distances from each other, the resisting pressure 
at a surface point of the body being proportional to the square of the normal 
component of its velocity, but without explaining how he obtained it. He didn't 
write, however, "I have a great proof, but no space for it in the margins of 
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this book" . A proof "from the Book" was waiting for the Pontryagin maximum 
principle. 

When one writes the resistance force TZ associated to Newton's problem, 



one obtains an integral functional of the type of those studied throughout the 
history of the calculus of variations. However, due to the restrictions on the 
derivatives of admissible trajectories, x{t) > 0, no satisfactory theory is avail- 
able within the calculus of variations framework (see 1, ,25 , 26 j. As first noticed 
by Legendre in 1788 (see [2] and references therein), without such restrictions 
on the derivatives the problem has no solution (the infimum is zero), since 
one can obtain arbitrarily small values for the integral resistance TZ [x(-)] by 
choosing a zig-zag function x(-) wildly oscillating, with large derivatives in ab- 
solute value. To make the problem physically consistent one must take into 
account the monotonicity of the profile, and this means, as was first remarked 
by V.M. Tikhomirov (cf. QJEDj that Newton's problem belongs to optimal 
control: 



Most part of the literature wrongly assume Newton's problem to be "one of the 
first applications of the calculus of variations" but, in spite of this, the same 
literature correctly asserts the birth of the calculus of variations: 1697, the 
publication date of the solution to the brachystochrone problem, and not 1686, 
the publication date of the solution to Newton's problem of minimal resistance. 

In 1997 H. J. Sussmann and J. C. Willems, in the beautiful paper de- 
fended the polemic thesis that the brachystochrone date 1697 marks not only 
the birth of the calculus of variations but also the birth of optimal control. The 
truth seems to be deeper: optimal control was born in 1686, before the calculus 
of variations, with Newton's problem of minimal resistance. The restriction on 
the control it(-), which appear in Newton's problem is a common ingredient 
of the optimal control problems. Such constraints appears naturally in practical 
engineering control problems, and are treated with the Pontryagin maximum 
principle - the central result of optimal control theory, first conjectured by 
L. S. Pontryagin, and then proved, in the late 1950's, by him and his collabo- 
rators, V. G. Boltyanskii, R. V. Gamkrelidze, and E. F. Mishchenko 22 . In an 
optimal control problem, the control functions take values on a set which is, in 
general, not a vector space. This is precisely what happens in Newton's classi- 
cal problem and the reason why Newton's problem must be classified as an 
optimal control problem, and not as a problem of the calculus of variations. 

Newton's problem has been widely studied, and the literature about it is 
extensive. The main difficulty is that of existence |S]: the Lagrangian L(t,u) = 





x(t)=u(t), u(t)>0, 
a;(0)=0, x(T) = H. 



(1) 
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jT^p associated to Newton's problem (Q is neither coercive nor convex, and 
Tonelli's direct method (see ^U]) fails. In order to prove existence, several 
different classes of admissible functions have been proposed. The question is 
now usually treated with the help of relaxation techniques (see jSJ), although 
direct arguments are also possible (see ^3 El)- As we shall prove (©, for 
the Newton-type problems, the existence of a minimizer follows easily from the 
Pontryagin maximum principle: one can show that the Pontryagin extremals 
are, for such problems, absolute minimizers (cf. Theorem EJ. 

Several extensions of Newton's problem have been considered in recent years. 
This revival of interest in Newton's problem, and in the study of many varia- 
tions around it, has been motivated by the paper j^j of G. Buttazzo and B. Ka- 
wohl. Recent results on Newton-type problems include: bodies without rota- 
tional symmetry (nonsymmetric cases) unbounded body (resistance 
per unit area) with one- impact assumption |12| : bodies with rotational symme- 
try and one- impact assumption, but not convex [U]; friction between particles 
and body (non-elastic collisions) ^3]; bodies with prescribed volume [3J; multi- 
ple collisions allowed ^JE]; unbounded body and multiple collisions allowed 
|20| . More recently, Newton-type problems have been related with problems of 
mass transportation |18l I19j . For a good survey on mass optimization problems 
and open problems, we refer the reader to [Sj. 

Here we consider convex <i-dimensional bodies of revolution with Height H 
and radius of maximal cross section T, and treat them using an optimal control 
approach. We will not be restricted to two-dimensional or three-dimensional 
bodies, considering bodies of arbitrary dimension d > 2. We also introduce a 
different point of view. For us the body does not move, and the particles are the 
ones who move. The body is situated in a flux of infinitesimal particles, the flux 
being invariant with respect to translations and rotations around the symmetry 
axis of the body. This new point of view is, in our opinion, physically more re- 
alistic. Newton has considered the particles with no temperature (not moving) . 
When the particles have temperature, they move, and the flux of particles is 
not necessarily falling vertically downwards the body, as considered by Newton. 
We will be considering new interesting situations with a non-parallel flux of par- 
ticles. We obtain complete solution to this class of Newton-type problems, by 
showing that, under some physically relevant assumptions on the Lagrangian, 
a control is an absolute minimizing control for the problem if, and only if, it 
is a Pontryagin extremal control. Thus, for the Newton-type problems we are 
dealing with, Pontryagin maximum principle holds not only as a necessary op- 
timality condition, but also as a sufficient condition. As very special situations, 
one obtains the solution found by Newton himself f H4.2.2fl . and solutions to 
Newton's problem in higher-dimensions ( ^4.2.3(1 . 
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2 Optimal Control 



The optimal control problem in Lagrange form consists in the minimization of 
an integral functional 

J[x(-),u(-)} = [ L{t,x{t),u(t))dt (2) 
Jo 

among all the solutions of a differential equation 

x'(t) = tp(t,x(t),u(t)) , te[0,T] (3) 

subject to the boundary conditions 

x(0) = a, x(T)=/3. (4) 

The Lagrangian L and the velocity function ip are defined on [a, b] x R n x f2, where 
fl C R r is called the control set. The main difference between the problems of 
optimal control and those of the calculus of variations, is that £1 is in general not 
an open set. In the case ip(t, x, u) = u, and Q = R", one gets the fundamental 
problem of the calculus of variations. For the Newton problem, we have n = 
r = 1, fi = Rq", <p(t,x,u) — u, a — 0, [3 > 0, and L(t,x,u) = 1+ * Mii . Typically, 
L(t,x,u) and tp(t,x,u) are continuous with respect to all arguments and have 
continuous derivatives with respect to x; the admissible processes (x(-),u(-)) 
are formed by absolutely continuous state trajectories x(-) and measurable and 
bounded controls u(-), taking values on the control set fl and satisfying ©-0- 
The Pontryagin maximum principle is a first-order necessary optimality con- 
dition, which provides a generalization of the classical Euler-Lagrange equations 
and Weierstrass condition, to problems in which upper and/or lower bounds are 
imposed on the control variables. 

Theorem 1 (Pontryagin maximum principle). Let (x(-),u(-)) be a mini- 
mizer of the optimal control problem. Then there exists a pair (tpo,ip(-)), where 
%j)Q < is a constant and ip(-) an n-vector absolutely continuous function with 
domain [0, T], not all zero, such that the following holds true for almost all t on 
the interval [0, T] : 

(i) the Hamiltonian system 
x'{t) = §£(t,a;(f),u(*),Vo,lKt)) , 

(ii) the maximality condition 

H (t, x(t), u(t), ipo, ^(*)) = maxH (t, x(t),u, ip ,i/j(t)) ; (5) 

where the Hamiltonian H. is defined by 

Tt(t, x, u, tpo, ip) = ipoL(t, x, u) + ip ■ ip(t, x, u) . 
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The first equation in the Hamiltonian system is just the control equation 
The second equation is known as the adjoint system. 

Definition 2. A quadruple (x(-), u(-), i/jq, ip(-)) satisfying the Hamiltonian sys- 
tem and the maximality condition is called a Pontryagin extremal. The control 
u(-) is said to be an extremal control. The extremals are said to be abnormal 
when ?/>o = and normal otherwise. 

Remark 3. If (x(-),u(-),ipo,il>(-)) is a Pontryagin extremal, then, for any 7 > 0, 
(x(-), u(-), 7^0) T0(")) ^ s a l so a Pontryagin extremal. From this simple obser- 
vation one can consider, without any loss of generality, that tpo = — 1 in the 
normal case. 

Remark 4. The fact that Theorem Q] asserts the existence of Hamiltonian mul- 
tipliers tpo ari d ip(') n ot vanishing simultaneously is of primordial importance: 
without this condition, all admissible pairs (&(•), would be Pontryagin ex- 
tremals. 

In some situations, it may happen that functions L and/or tp depend upon 
some parameters w S W C R fc . In this case, given a control u(-), the corre- 
sponding state trajectory x(-) and the cost functional J will in general depend 
on the choice of the parameters w. The problem in then to choose the param- 
eters w in W for which there exists an admissible pair (x(-),u(-)) such that 
J [x (•),■&(•), w] < J [x(-), u(-), w] for all w € W and corresponding admissible 
pairs (x(-),u(-)). The parameter problem can be reformulated in the format 
©-® by considering w as a state variable with dynamics w'(t) = and initial 
condition w(0) € W. 

3 Optimal Control of Newton- Type Problems 

The standard method to solve a problem in optimal control proceeds by first 
proving that a solution to the problem exists; then assuring the applicability 
of the Pontryagin maximum principle; and, finally, identifying the Pontryagin 
extremals (the candidates) . Further elimination, if necessary, identifies the mini- 
mizer or minimizers of the problem. It is not easy to prove existence for Newton's 
problem with the classical arguments, because the Lagrangian L(t, u) = is 
not coercive and it is not convex with respect to u for u > 0. Here we will make 
use of a different approach. We will show, by a simple and direct argument, 
that for Newton-type problems ©-0 the Pontryagin extremals are absolute 
minimizers. This means that, in order to solve a Newton-type problem, it is 
enough to identify the Pontryagin extremals (cf. Theorem [BJ. 

We begin to show that there are no abnormal extremals for a Newton-type 
problem. 

Proposition 5. Let L(t,u) be a continuous function satisfying the following 
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conditions: 

L(t,u)>£>0 V(t,u) €]0,T] x M+ , (6) 
lim L(t,u)=£ Vt€ [0,T]. 

u — > + oo 

Tften aZZ Pontryagin extremals (x(-),u(-),ipo,ip(-)) of the Newton-type problem 

J[u(-)]= [ L(t,u(t))dt — >min, 
Jo 

x'(t)=u(t), u(t)>0, (7) 
x(0) =0,x(T) = (5 with [3>0, 

are normal extremals (ipo = — 1 ) with tp a negative constant (ip(t) = —A, A > OJ. 
Proof. As far as the Hamiltonian does not depend on x, 
7i (t, u, ip , ip) — ipoL(t 7 u) + tpu , 

we conclude from the adjoint system that tp(t) = c, with c a constant. If c 
is equal to zero, then i/jq < (they are not allowed to be both zero) and the 
maximality condition JSJ simplifies to 

tpoL(t,u(t)) = max {ipoL(t, u)} . 

Under the hypotheses © the maximum is not achieved (u — > +oo) and we 
conclude that c ^ 0. Similarly, for c > the maximum 

max {ipoL(t, u) + cu} 

does not exist and one concludes that c < 0. It remains to prove that tpo is 
different from zero. Indeed, if ipo — 0, © reads 

cu(t) = max {cm} , 

U>0 

and follows that u(t) = and x(t) = constant. This is not a possibility since 
(3 > 0. The proof is complete. □ 

Theorem reduces the procedure of solving a Newton-type problem to the 
computation of Pontryagin extremals. 

Theorem 6. The control u(-) is an absolute minimizing control for the Newton- 
type problem ©-(0, i.e., J[u(-)] < J[u(-)] for allu(-) € ([0,T],M+), if and 
only if, it is an extremal control. 

Proof. Theorem is a direct consequence of the maximality condition. From 
Proposition [5] one can write JSJ as 

-L(t, u(t)) - Xu(t) > -L(t, u(t)) - Xu(t) (8) 
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for all admissible controls u(-) and for almost all t £ [0, T\. Having in mind that 
all admissible processes (x(-),u(-)) of the problem J7J) satisfy 



The required optimal solutions of the Newton-type problem are ex- 

actly the Pontryagin extremals. This means, essentially, that we have reduced 
a dynamic optimization problem (a minimization problem in the space of func- 
tions) to the static optimization problem given by the maximality condition. 

Corollary 7. Finding the solutions for the Newton-type problem amounts 
to find the minimum of the function h{u) = L(t, u) + Xu, t £ [0, T], A > 0, for 
u>0. 

4 An Application 

Consider a d-dimensional body of revolution 



where d > 2, £ = . . . , £d-i), $ is a non- negative function defined on 
[0, if]. Denote by T the radius of maximal cross section of the body, T = 
maxo<^<H 3>(C)- Let us assume that the body is convex, then the function $ is 
concave, and there exists c £ [0, if] such that $(C) is monotone increasing as 
C < c, and monotone decreasing as £ > c. 

We suppose that the body is unmovable and is situated in a flux of infinites- 
imal particles. The flux is invariant with respect to translations and rotations 
around the £o-a-xis, which is the symmetry axis of the body. So, the specific 
pressure of the flux on an infinitesimal element of the body surface depends 
only on the value of <&' at that element. It is convenient to consider, instead of 
<£>, two functions that are generalized inverses of <&; denote them by X-(t) and 
H—x+(t), They are defined in the following way: x_ (t) ~0&st£ [0, $(0)], and 
X-(t) is inverse to the strictly monotone increasing branch of $ as t £ [^(O), T]; 
x + (t) = as t £ [0, <f>(if)], and if — x + (t) is inverse to the strictly monotone 
decreasing branch of $ as t £ [$(if), T]. The obtained functions x^. and x + 
are convex, continuous, and monotone increasing, besides x_(0) = x + (0) = 0, 
X-(T) < c, x+(T) < if — c. In such a representation, the specific pressure is 
a function of x' + or of x'_ , if the point belongs to the front or to the rear part 
of surface, respectively; we denote the corresponding functions by p+{ ) and 




it is enough to integrate ||SJ| to obtain the desirable conclusion: 




□ 



{(£o, ■ Co G [0, H], |f| < $(£ )} C R d , 



*-(•)• 
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The pressure on an element d d 1 s of the front part of surface is dp + = 
p+(x'_i_)d s. The projection of the pressure vector to the £o-& x is equals dpo = 

dp+ 1 yjl + x'f , and the projection of the surface clement to R^ -1 ^ 1 has area 

d^ 1 ^ = d d ~ 1 s/ ^Jl + x'f. Thus, the £o -projection of pressure corresponding 

to the element is dpo = p+(x' + )d d ~ 1 ^. Passing to polar coordinates and 

integrating over the ball {|£| < T}, one obtains the resistance TZ+ of the front 
part of body to the flux: 

n + [ x+ (-)\ = [ T P+ {x' + (t))dt d -\ 

Jo 

here Vd-i stands for the volume of (d — l)-dimensional unit ball. Similarly, 
the resistance of the rear part of body to the flux (which is positive) equals 
— TZ-[x _(•)], where 

H-[x-(-)] = v d - 1 f P -{x'_{t))dt d -\ 

Jo 

So, the resistance of body to the flux is TZ[x+(-), X-(-)] = 1Z+ [x+(-)] +1Z- [x -(■)]. 

It is required to minimize TZ[x+(-), X-(-)] over all pairs X-(-)) of 

convex monotone increasing functions defined on [0, T] , provided x± take values 
in [0, (3±], where /?_ = c, (3 + = H — c, T and H are fixed, and c varies between 
and H. 

We are acting as follows. First wc fix the sign "+" or "— ", minimize 1Z± 
over monotone increasing functions x : [0, T] [0, [3±], with (3± fixed, and 
verify that among all the solutions, the convex one is unique; denote it by x^± . 

Then we minimize the sum TZ+[x^ (■)] + lZ-[x_ (•)] over all positive (i + and 
(3- such that /?++/?" = H. 

4.1 Solving the problem in general case 

In what follows, we assume that the functions p + and p- satisfy the following 
conditions: 

(i) PiSC^O, +oo); 

(ii) there exist lim„^ +oc p±(u); 

(iii) p' ± (0) = lim„^ +oc p' ± (u) = 0; 

(iv) for some u± > 0, p' ± is strictly monotone decreasing on [0, u±], and 
strictly monotone increasing on [u±, +oo). 

For simplicity, we further put Vd-i = 1. Let us fix the sign "+" or "— ", and 
introduce shorthand notations 

p± =p, P± = 13, n± = n, xi± = x p. 

Proposition 8. There exists a unique solution u° of the problem 

p(0)-p(u) 

► max, 

u 
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besides u° > u. 

Proof. Denote q(u) = p(0) — p(u) and B = sup u>0 (q(u) / u) . From (i)-(iv) it 
follows that the function q(u)/u, it > is continuous, positive, and satisfies 
the relations lim. u ->o+(q(u)/u) = lim u ^ +00 (q(u) /u) = 0, hence < B < +oo, 
and there exists a value u° > such that q(u°)/u° = B. Obviously, at u = u° 
one has (q(u)/u)' — 0, hence q'(u°) — q(u°)/u°. At some 9 £ (0, 1) one has 
q(u°)/u° = q'(0u°), hence q'(u°) = q'(8u°). This implies that q' is not strictly 
monotone on [0, u ]; thus, by virtue of (iv), vP > u. 

It remains to prove that the value u , solving the equation q(u)/u = B, 
is unique. Suppose that q(u°)/u° = q(u )/u — B, u° < u . Then q{u°) = 
u°q'(u°), qiu 1 ) — u 1 q'(u r ). At some u £ (u°, u 1 ), one has q(v}) — q(u°) — 
q'(u) (u 1 — u°): this implies that 

q'{u) {u 1 - u°) = u 1 q'iu 1 ) - u° q'(u°), 

hence 

u° (q'iu ) - q'{u)) + u 1 [q'{u) - q' [u 1 )) = 0. (9) 

One has u° > u, hence q' is strictly monotone decreasing as u > u , so both 
terms in are positive. The obtained contradiction proves the proposition. □ 

Let us denote 

p(0)-p(u) , 

B = = — p (u). 

u 

Proposition 9. (a) As \t 2 ~ d > B, the unique solution of the problem 

t d ~ 2 p(u) + Xu^ min; (10) 

is u = 0. 

(b) As At 2 ~ d — B, there are two solutions: u = and u = u° . 

(c) As \ t 2 ~ d < B, the solution u is unique, besides u > u° , and p'(u) = 

-xt 2 - d . 

Proof, (a) and (b) are obvious; let us prove (c). Denote A := \t 2 - d . By 
definition of B, for < u < u° one has 

P(0)-P(") ^ o P(0)-P(u°) 
<B = 5 > 

p(u°) + Bu° = p(0) < p(u) + Bu , 

hence 

p{u) - p{u°) >B{u°-u)>\ (u° - u), 

and thus, 

p{u) + Xu> p(u°) + Xu°. 
On the other hand, one has B = —p'(u°), therefore 

p'(u°) + A < 0; 
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moreover the function p(u) + Xu is convex on [u°, +00) and tends to infinity as 
u — > +00. All this implies that the solution u of IjlOl) is unique, satisfies the 
equation p'(u) + X = 0, and u > u°. □ 

From Corollary we know that if x^(-) is a solution of the minimization 
problem 72.[x(-)] — > min, x : [0, T] — > [0, /3], then for some A, the values 
u = x 13 (t), t E [0, T] satisfy the equation (|10|1 . According to propositions 
121 and [21 one should distinguish between three cases: (a) if Xt 2 ~ d > B, then 
u = 0; (b) if Xt 2 - d = B, then u = or w = u°; (c) if A< 2 " d < S, then it > u°, 
and p'(u) = -Xt 2 - d . 

Consider two different cases: d — 2 (two-dimensional problem) and d > 3 
(the problem in three or more dimensions). 



4.1.1 Two-dimensional problem (d = 2) 

If A > B, the unique solution of (|10|l is u = 0, hence x 13 = 0. This implies that 
,9 = 0. 

If A = B, there are two solutions: u = and u = u , therefore any absolutely 
continuous function x(-), x(0) — 0, x(T) = [3, such that x'(t) takes the values 
and u°, minimizes 7Z. A convex solution x@ has monotone increasing derivative, 
hence for some to G [0, T], x 13 (t) = as t G [0, to], and x 13 (t) = u° as 
t e [to, T). Thus, 

(3 (f] _ f aste [0, t ] 

1 1 j ~ \ u°(t-t ) aste [t , T]. liiJ 

Taking into account that x@{T) = (3, one concludes that (3/T < u° and to = 
T-0/u°. 

If A < B, there is a unique solution u, hence ar (t) = ut, and /3/T = u > u°. 



Summarizing, one gets that 

(i) as (3/T < u°, the convex solution x"(t) is given by (|llf> : 

(ii) as [3/T > u°, xP(t) = /3t/T. 

As (3/T < u°, one has 

R[xP(-)]= [°p(Q)dt+ f p(u°)dt, 

JO J t 

and taking into account that t = T — (3/u° and (p(0) — p(u°))/u° = B, one gets 

ft[x"(-)] =Tp(0)-(3B. 
As /3/T > u°, one has U[x^(-)] = Tp((3/T). Introduce the function 



p(u) 



p(0) - Bu, if < u < u°, 
p(u), if u > u°, 
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p'±(u) = < , A. 



then 

= Tp(/J/T). 

Thus, the minimization problem lZ+{x^ (■)] + 7Z- [x^_~ (•)] — > min is reduced to 
the problem 

Ph(z) — p+(z) + p- (h — z) — » min, < z < /i, (12) 
where h = H/T. 

The introduced functions p±(u) are continuously differentiable on [0, +oo), 
and 

-B± if < u < u° ± , 

4- 

Using that it^_ > u±, one concludes that p'±(u) is monotone increasing, hence 
p' h (z), < z < h is also monotone increasing. 

From now and until the end of subsection l4.ll we shall assume that p' + (u) < 
p'_(u), u > 0, hence B + > £>_. Denote by it* a positive value such that 
p'Au*) = —B-. This value is unique, and it* > iq_. Consider four cases: 
1) < h< u+; 2) u\ <h< u*; 3) it* < h < it* + and 4) ft > it* + u° . 

In the cases 1) and 2) one has 

Ph(z) <Ph( h ) =P'+( h ) + B - ^ aS < Z < ^ 

hence the minimum of (|12[l is achieved at z = h. Therefore, the optimal value 
of j3- is zero, so x_ = 0. 

In the case 1) one has (3 + /T = h < i/! , hence x+ + (•) is given by (jllfl . with 
= T (1 — h/vP + ). So, the optimal body is a trapezium. 

In the case 2) one has x^ (t) = /it, hence the optimal body is an isosceles 
triangle. 

In the cases 3) and 4), one has p' + (h) > — £?_ > — B + = p' + (u+), hence 
h> vP + . Further, one has 

p' h (h)=p' + (h)-B- >0; 

on the other hand, 

p' h (u° + ) = - p'_{h ~ u° + ) < -B + + B-<0. 

It follows that the minimum of ph is achieved at an interior point of [u+ , h] , so 
the optimal value of /3_ satisfies the relation u°^_ < (3+/T < h, and x^~ (t) = 
t/3+/T. 

In the case 3), denoting h — max{0, h — u°_}, one has h < it*, hence 
p' h Ch)=p / + (h)-p'-(h-h)<p' + Ch) + B- <0, 

hence the minimum of ph is reached at an interior point of [h, h], thus < 

0-/T <h-h<u_, and 

b (0 ift€[0, T-/3_/u°] 

- u°_{t-T + l3-/u -) if t e [T-P-/u°_, T]. 
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The optimal body here is the union of a triangle and a trapezium turned over. 

In the case 4), one has p' h (h — u°_) = p' + (h — «_) + £?_ > 0, hence the 
minimum of ph is reached at a point of [11+, h ~ vP_). Thus, (3-/T > u°_, and 

x^T (t) = t/3-/T. The optimal body is a union of two isosceles triangles with 
common base. 



4.1.2 Problem in three or more dimensions (d > 3) 

Here we additionally assume that p± € C 2 [0, +00) and > as u > v? + . 

Denote lj = l/(d— 2) and to = (X/B)^. As < t < to, the unique solution 
of UlUfl is u — 0, hence ar(i) = 0. As to < t < T, the solution u satisfies the 
relation 

t d - 2 p'{u) + X = 0, 

and u > u. 

If T < to, one has x 13 = and j3 = 0. Let, now, to < T; using that 
p' is negative, continuous, and strictly monotone increasing on [u°, +00), one 
concludes that x 13 (t) = u is also continuous, and is strictly monotone increasing 
on [to, T)} from x 13 {to+) — u° to the value U defined from the relation 

T d - 2 p'(U) + A = 0, U>u°. 

Thus, x {-) is convex; moreover, as to < t < T, x 13 can be represented as a 
function of u G [vP , U\. Using that x 13 ' ' (t) = u and that 



(13) 



one gets 



hence 



dx 13 dx 13 dt , , , d ( 1 
-— = uA" — 



du dt du du \ \p'(u)\ u 

x f3 = \" / vd 



\p>(y)\> 



using that = B, one obtains 

u u° f u dv 



x 



f} = X u 



In particular, substituting u = U, one has 

A" ( U f U dV \-3 (U) 

[WW B» J u o \p>(y)\«) P - 1 ' 

Introduce the function 

, . r dv 

9{u) 



WW 
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Using that |p'(f)| = B as < v < u°, and p' iy) = p'{v) as v > u°, one gets 

, 7r , u° A 1 ' di/ 

9{u) = b^ + L RRP 

and using that 

T= A 



b'([/)K 
from l|14[l one gets 

| = f/-b'(c/)r 5 (f/). as) 

The minimal resistance equals 

,T ,T 

K[ X P(.)]= piuitydt*- 1 =p(Q)t d Q - 1 + p{u(t))dt d -\ 

JO J t 

Using that u(T) = £/, u(i ) = u°, |p'(u°)| = B, p(0) - p(tt ) = Bu", t Q = 
A"/-B w , and also the formula (fH^ . one finds 



0\ /■{/ 



dp(u) 



This implies 



^P=p(C7) + |p'(C7)|i+^(C/). (16) 



Denote E/ + = z + , [/_ = z_. Using i|15|) and (|16[) . one comes to the following 
problem of conditional minimum 

r(z-,z+) := 

p + (z + ) + |p' + (z+)| 1+ - g + (z+) + p_(*_) + -> min, 

under the conditions 

Z--\ P , _(z-)r 9 -(z-)+z + -\p'_(z + )\»g + (z + ) = h, z_> u °, z+> u °. (17) 

From l|17|l. taking into account that \p'±(z±)\ u> g'±(z±) = 1, one obtains that z+ 
is a differentiable function of z_, and 

= b / _(z-)|"- 1 p^(z-)g-(z-) 
d*- |pU^)I w_1 *»+(«+) S+(*f)' 
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Now, 

— r(z-,z+{z-)) = 

= - 1 ±.{l + u)\p' + {z + Yp'l{z + )g + {z + ) (l + u,)\p'_(z-Wp'L(z-)g-(z-) 

= (i+w) ij/_(z_)r-v^(2-)s-(«-) • -*/+(*+)). 

Note that z+ is a monotone decreasing function of z~, hence the function 
p'_{z-) — p',(z+(z—)) is monotone increasing as z_ > u , > u . 

Recall that u* is the value satisfying = — £?_. Denote 

/i* = w* — g+(w»). 

Consider two cases. 

1) /i < ft,*. One has z + (u°) < u„, hence f/_(tt°) — p',(z + (w )) > — £>_ — 
p'+(u*) = 0. It follows that as z_ > u , p'_(z~) — p' + (z + (z_)) > 0, so the 
minimum of r(z_, z+)is attained at z_ = 0. 

2) h > K. One has z + (u°) > u*, hence p'_(u°) - p' + (z + (u )) < 0. On the 
other hand, as z_ = z + (z_), one has p'_{zJ) — p'Az-) < 0, hence at some 
Z- € (it , p'_{z^) = p' + (z + (z_)), and so, the minimum of resistance is 
attained. 

4.2 Examples 

We have given in fl4.ll complete description of the solutions to the formulated 
Newton-type problem. We now consider, for illustration purposes, various par- 
ticular cases of the problem. All the calculations can easily be done with the 
help of a computer algebra system. We have used Maple to implement a proce- 
dure which, given functions p+{-) and P-(-) and the values for T and H, gives 
the optimal shape for the respective problem. 

4.2.1 Non-parallel flux of particles 

Let us consider the two-dimensional case (d = 2). As proved in H4.1.1I there exist 
four possible cases. To illustrate this we choose, as an example, the pressure 
of the front part of the surface to be p+ — jt^t + 0-5; the pressure on the 
rear part given by p- — fr^p — 0.5; the radius T of the maximal cross section 
of the body to be two (T = 2); and then we choose different values for the 
height H of the body. Applying the formulas given in i!4.1.1l one obtains that 
for H = 1 the solution is a trapezium (Fig. QJ; for H = 2 a triangle (Fig. 0); 
for H = 4 the union of a triangle and a trapezium turned over (Fig.^J); and for 
H = 6 the union of two triangles with common base (Fig. E}. We remark that 
in Newton's problem one has p+ = 1 _^ u2 and p- = (parallel flux), and only 
the first two situations occur: solution to Newton's two-dimensional problem is 
either a trapezium or a triangle. 
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Figure 1: H = 1 Figure 2: H = 2 




-2-1 12 -2-1 12 

t t 



Figure 3: H = 4 Figure 4: = 6 

Solutions of the two-dimensional Newton-type problem with p + = yq^u + 0.5, 
P- = TTm 2 ' ~ (non-parallel flux of particles), T — 2, and different values for 
the height H of the body. 
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The two-dimensional problem under a non-parallel flux of particles with 
density of distribution over velocities circular gaussian, with biased mean, is 
studied in |2"T] . 



4.2.2 Newton's classical problem 

We now obtain the well-known Newton's solution. For that we fix d = 3, 
p+(u) = 1/(1 + it 2 ), and p-{u) = 0. Applying the method described in H4.1.2I 
after some algebra one obtains u + = 1/v3j u ° = 1, B + = 1/2, (3 = H, and 
the optimal solution x(t) is given in parametric form by 

A /3u 4 2 , 7 

x = — hw — In u 

2 V 4 4 

t = - f u 3 + 2u + - ) , 1 < u < U . 
2 V uj 

all in agreement with classical formulas. Expressing the formulas with respect 
to U and T one obtains: 

2TU t ATU _ TU (-7 + 4U 2 + 3U 4 - 41n(U)) 



"(1 + C/ 2 ) 2 ' ~(l + t/ 2 ) 2 ' ' 4 (1 + U 2 ) 

TU(l + u 2 ) 2 TU (-7 + 4it 2 + 3u 4 - 41n(w)) 

u(l + U 2 ) 2 ' 4(l + [/ 2 ) 2 ' 

_ T 2 (17J7 2 + 2 + 10C/ 4 + 3J7 6 + 4 ln([/)[/ 2 ) 
4(l + (7 2 ) 4 

In this case 7£_ = 0. 

4.2.3 Newton's problem in higher dimensions 

Our approach to Newton's problem is valid for an arbitrary d > 2. For example, 
for d — 4 (problem in dimension four) one gets: 



2T 2 U 



^ - T. 77^2 > *0 — \ / — — ttt > - 



4T2[/ T (-5 17 + 3 U 3 + 2VU 



(1 + U 2 ) 2 ' ' \l{l + U 2 ) 21 r 5(l + U 2 ) 

t = T [U(l + u 2 \ ^_TVU(-5^I+3u 5 / 2 + 2) 



u \1 + U 2 J ' 5(1 + U 2 ) 

T 2 (l + 3U 2 ) 



2 (1 + U 2 ) 2 
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